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Abstract 

We show that when the magnetizations on opposite sides of a junction are not collinear, the 
magnetization current is not conserved as the junction is crossed. Thus the usual treatment of this 
problem needs to be modified. We argue that this is due to an implicit assumption of an external 
torque that is required to clamp the magnetization in place. The physical consequence of this is 
explored. 

PACS numbers: PACS numbers:73.40.-c,71.70.Ej,75.25.+z 



There is much recent interest in spin polarized transport in heterostructures such as 
tunnel junctions and metallic multilayers. A general approach to these types of problems 
is to solve for the steady state solutions. To incorporate the transmission and reflection of 
the injecting current at the interface, it is necessary to solve the boundary value problem. 
Central to treating the transmission and the reflection is how the charge current J and 
magnetization current Jm go from one side of the heterostructure to the next. For a simple 
heterostructure with no magnetic impurities both types of currents are usually assumed to 
be conserved as they go across the barrier. An equation that is usually employed in the 
treatment of this phenomena is of the formQ] 

Jm = ^(as^Ls - PsVRs) (1) 
s 

where a, (3 are coefficients, Pl(r)s are the chemical potentials of the electrons of spin channel 
s on the left (right) of the junction. In this equation, it is assumed that Jm is unique 
because it is conserved. Recently, it was discovered that when the magnetizations on opposite 
sides of a metallic barrier structure are noncollinear, a current can induce torques on the 
magnetizations . This motivated us to examine the boundary condition across a 

heterostructure more closely. We consider the strong barrier limit corresponding to tunnel 
junctions and the weak barrier limit with zero interface resistance, corresponding to metallic 
multilayers. We find that when the magnetizations on opposite sides of the barrier are 
noncollinear, the magnetization current is not conserved as it goes across the barrier. Thus 
eq. (1) cannot always be used. An examination of the physics suggests that implicit in the 
assumption of the noncollinear magnetizations is a torque to maintain the orientation of the 
magnetization. It is this torque which changes the magnetization current. We now discuss 
our results in detail. 

We first remind the reader the meaning of the direction of magnetization. The general 
quantum statistical properties of a system is described by the density matrix p. The general 
state of magnetization is therefore described by the 2x2 polarization matrix |Mo|p where 
\M \ is the magnitude of the magnetization. As is well known[6] there is always a direction 
M so that p = (1 + a • M) where a are the Pauli matrices (a x , a y , a z ). The density matrix 
is diagonal along M. We call this the direction of magnetization. Along this direction, 
the off-diagonal components of the density matrix are zero and there is no phase coherence 
between the spin up and the spin down states. 



We next evaluate the magnetization current. For simplicity we take a one dimensional 
geometry where the currents are uniform in the x-z plane and move in the y direction. 
The evaluation of the magnetization current is simplest when we pick a coordinate system 
along the direction of magnetization. The magnetization current then becomes the difference 
between the spin up and the spin down current. We first consider the strong barrier limit 
of spin polarized tunnelling. Our starting point is the well accepted tunnelling Hamiltonian 
Ht = 'YliT ss ic[ ks cn qs i + c.c. where s, s' are with respect to the magnetization directions on 
the left (labelled by the subscript L) and on the right (labelled by the subscript R). k, q 
are other labels for the electronic states. We are interested in the case when the directions 
of magnetization on opposite sides of the junction are not the same. The spins states \s) 
along one magnetization direction are related to those {\s')) along the other direction by a 
rotation matrix R SjS /. From a single particle quantum mechanical calculation we obtain an 
estimate of the tunnelling matrix element 

rd 

T ss , = R ss , exp[ / dy{h\U - \i La + ^ - Vy)/2m) - 5 }. (2) 
J o 

where U is the tunnelling barrier, V is the external voltage, m is the electronic mass. In 
linear response, the term — \i Ls + \i Rs i — Vy is absent. 

Stimulated by the experimental success in superconductivity tunnelling there has been 
much studies done with this type of tunnelling Hamiltonian. We shall follow the work by 
Cohen, Falicov and Phillips!?], as is explained in detail in Mahan^. The magnetization 
current on the left can be evaluated from the time rate of change of the left magnetization 
due to the tunnelling Hamiltonian: Since J s = —dN s /dt, we get Jj^ = —d(N+ — Nt)/dt = 
—dMi/dt (we have used a unit so that the Bohr magneton is 1) with a similar expression for 
Jpji- The time derivative of the particle number of spin s can be evaluated as its commutator 
with the tunnelling Hamiltonian. The details of this for each spin component is essentially 
the same as in previous calculations and thus will not be repeated here. We find thatj^] 

J\ L = 47re 22 s\Tss'\ 2 {n L (^ Ls ) - n R (n Rs <)), (3) 

ss' 

and 

J§ = Ane 2j s'\Tss'\ 2 {n L (^ Ls ) - n R {n Rs <)). (4) 

where hl s i^Rs 1 ) is the chemical potential for the spin s (s') electrons on the left (right) side. 
ti(ijl) is the density of particles at energy \i. 



When the directions of magnetization on both sides of the junction are collinear, T is di- 
agonal and we find that J M = 4ne J2 S s \Tss\ 2 (riL(fiLs) — Wr(jUrs)) = Jm- The magnetization 
current is conserved. When the magnetization directions are not collinear, T is no longer 
diagonal. In the linear response regime, T SjS i = T S > )S is a symmetric matrix. After renaming 
the variables of summation, one has J M = 47re^ ss , s'\T 'ss'\ 2 '{n^ni,^) — n R (fj, Rs )). J M ^ J M 
because hl s ' ^ Hl s for s ^ s'. The latter comes about because there is a splittin g of the spin 
up and spin down chemical potential caused by the spin accumulation effect jloL ll l| and the 
interplay between the charge and magnetization dipole layers [121]. 

J M and J M correspond to magnetizations along different directions. In the above we 
have shown that the magnitude of these two currents are not the same. A more appropriate 
condition would be requiring the magnetization currents for magnetizations along the same 
direction be equal: J M = J M . This requirement obviously cannot be satisfied because for 
example, on the right hand side it would require a magnetization current with a magneti- 
zation component that is perpendicular to the magnetization on the same side. From what 
we just derive, this is zero. The reason for the nonconservation of the magnetization current 
is easy to understand. Consider a wave function \sL) of an electron of spin s on the left. 
With respect to the magnetization direction on the right, this wavefunction can be written 
as J2 S ' ^ss'\s'R). The corresponding density matrix for a state like this is given by 



Po 



K ++ K + 



with non-zero off-diagonal elements. But by our definition of the magnetization direction, 
we have assumed that after this electron has moved from the left to the right, the phase 
coherence between the up and down states are lost and the off-diagonal components become 
zero. This happens because we have implicitly assumed that there is some external torque 
that clamps the magnetization in certain directions. This torque can come from the crys- 
talline anisotropy energy of the system, for example. It is this implicit torque that causes the 
magnetization current not to be conserved. In addition to this torque to = {Jm ~ ^mV^-B 
(fiB is the Bohr magneton, a is atomic spacing) there may be other contributions which 
we shall not discuss in this paper. The reasoning in this section also applies to metallic 
multilayers. We turn our attention to that case next. 

For metallic multilayers, eq. (3) and (4) are no longer valid. For simplicity we consider 
the other limit where the interface resistance is zero. For that case, the current is controlled 



by the voltage drop along the metal. We assume that the length of heterostructure is less 
than the spin diffusion length and apply a parallel circuit model to illustrate the essential 
physics. In that case, we get = crf^ E L ^ where crf^ is the conductance for spin 

s on the left (right). E is the electric field that drives the current. The spin current on 
opposite sides of the heterostructure are related by J^, = J2 S |R SjS /| 2 J^, where, again, R is 
the rotation matrix. There are two ways that one can see this. The current on the right for 
spin index s' is just the number of electrons with z angular momentum s' moving past per 
unit time per unit area. These electrons come from the left side. The factor R 2 takes into 
account the probability that the z angular momentum is s'. Formally, in terms of electron 
creation and destruction operator the current operator is just J R = (Yl vc t' c s') where v is 
the velocity. Now the creation and destruction operator can be expressed with the z axis 
along Ml as c s i = R s / jS c s , = R s / jS c+. Because R is real, it occurs in the expressions 
for both the creation and destruction operator. By our assumption of the magnetization 
direction, we have set such terms as < Yl c t c t > with s ^ t to be zero. Hence the final 
expressions quoted above is obtained. 

It is straightforward to show that E L( ~ RS} = a L ^E/(a L + a R ) where the total conductance 
a L(R) _ ^2 s (js^ R \ E = E L + E R measures the total voltage drop. Just as explained in the 
previous paragraph, because the magnetization points in different directions on opposite 
sides of the heterostructure, the magnetization currents are not equal to each other. The 
magnitudes of the two currents are equal to each other, however. 

In the presence of a finite interface resistance, in linear response the current is linearly 
proportional to the chemical potential. We thus expect the physics to be governed by 
equations of the form 

# = £r^-r£ty? (5) 

s 

# = £r£^-r**,x? (6) 

s 

where the r's are the effective interface resistances. We next explore the physical implication 
of our result. 

In typical experimental arrangements, only the directions of the magnetizations far away 
from the interface are specified. Because of the torque, these directions will be rotated as 
the interface is approached. A domain wall is created. We estimate this rotation here for a 
simple case. We assume that far from the interface the magnetzations are in the x-z plane 



along the easy axis of anisotropy at angles 9l, Or so that 9r — 6l = At the interface, 
the left (right) magnetization has rotated by an angle — 59 r (59 l). From minimizing the 
energy, the equation governing the angle 59 at a distance y away from the interface is 
given by Jid 2 59/dy 2 — Ki59 = where i = L,R; Jj ,Ki are the exchange and anisotropy 
constants. The effective field H e ff due to the exchange, Jid 2 59/dy 2 , is balanced out by 
that due to the anisotropy, K, b 59. The solution of this equation is 59(y) = 59iexp(—\y\/di) 
where the magnetic length d\ = \J Ji/Ki. At the interface, the effective field from the 
exchange becomes Jd59/dy while that from the anisotropy remain the same. The net torque 
t = -dE/d9 = M 59i(Ki - Ji/di) = M 59i[Ki - (J^) ' 5 ]. This creates the difference in the 
magnetization current discussed above. We thus get t = t . 59 t = AJ M / M a[Ki — ( Jji^j) - 5 ]. 
Thus we can relate 59 i to AJm- 

Experimental measurements are carried out under constant charge current conditions. It 
is necessary to relate AJm to the current J. To achieve this goal it is necessary to find 
solutions for the currents for each side of the junction and match them across the junction. 
In the treatment of spin polarized transport, we should use eq. (3) and (4) instead of eq. 
(1). These equations will be complemented by equations which relate the currents to the 
chemical potentials inside the material on either side of the junction. An example of such a 
equation is the modified Landau Gilbert equation with a source term equal to the divergence 
of the magnetization current: 

dM cl 5M 

-- 7 MxH + V-J m = 7 

at r 

where 7 is the gyromagnetic ratio, Jm is the spin current (tensor), r is the spin relaxation 
time, and H is the effective field describing the precession of the magnetic moments. The 
magnetization and charge currents are given by Fick's law as sums of terms proportional 
to the gradients of the charge (5n) and magnetization (5M) densities and the local internal 
electric field. They are 

J = —aW — f3DV5M ■ p — aiDVSn (8) 
j M = -/foWpo - DV5M - a 2 (3DX75np (9) 

where a,\ and a 2 are two phenomelogical parameters, (3 measures the asymmetry of the spin 
up and spin down current in the metal, a is the conductivity of metal, and D is the diffusion 
coefficient. p is the unit vector along the direction of magnetization. V is the sum of the 



electic potential V e describing the external electric field and the local electric (screening) 
potential W due to other electric charges determined self-consistetly 



W(v) = J dh'U{v-v')5n{v'). 



From these equations and the charge conservation 

„ . d5n , , 

v-J—a (10) 

we obtain for each side of the heterostructure the charge and magnetization dipole layers at 
the interface given mathematically by 

(11) 

i=i 

where the four length scales correspond to the screening length and the magnetic lengths for 
the longitudinal and left and the right rotating transverse magnetizations, d is the thickness 
of the barrier. The superscript R denotes the right hand side. A similar expression can be 
written for quantities on the left hand side. Because of the coupling between the charge and 
spin degrees of freedom, the length scales k are renormalized from their bare value. The 
coefficients 5riio are found to depend on 5~M i0 , while 5~M i0 are determined by matching the 
boundary condition at the interface. We illustrate this approach with a calculation of the 
steady state solutions in the low field linear response limit. After this boundary matching 
is carried out, we found for metallic multilayers in linear response a large change in the 
longitudinal component of the magnetization that is of the order of the voltage times the 
density of states at the Fermi surface (in units with /ib = 1)- The transverse magnetization 
change is much smaller. As a result, there is no additional torque generated along the 
ferromagnet. All the torque comes at the interface. 

When the material parameters on both sides are the same, we found 

eAJ M = /3Jcos($/2) sin$ 
M + 2 yT^fll + A 2 sin 2 (<l>/2)/(l - (3 2 )] + £ 2 4 sin 2 g/2 (12) 

1 - s- ■ 2 V /T^^c • c^in 2< l' 

where £ is the ratio of the conductance of the interface to that of the metal. A 1: A 2 are 
some dimensionless constants of the order of unity. 



In summary, we argue that the magnetization current is not conserved at the interface 
of a heterostructure when the magnetizations on opposite sides of the structure are not 
collinear. This is due to the implicit assumption of a torque which creates a domain wall 
at the interface. This is illustrated with a steady state solution of the of the problem in the 
limit of linear response. In that limit, all the torque comes at the interface, no additional 
torque is introduced in the metal. 
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